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Abstract 

In this paper, we construct an adiabatic invariant for a large 1-d 
lattice of particles, which is the so called Klein Gordon lattice. The 
time evolution of such a quantity is bounded by a stretched exponential 
as the perturbation parameters tend to zero. At variance with the 
results available in the literature, our result holds uniformly in the 
thermodynamic limit. The proof consists of two steps: first, one uses 
techniques of Hamiltonian perturbation theory to construct a formal 
adiabatic invariant; second, one uses probabilistic methods to show 
that, with large probability, the adiabatic invariant is approximately 
constant. As a corollary, we can give a bound from below to the 
relaxation time for the considered system, through estimates on the 
autocorrelation of the adiabatic invariant. 

1 Introduction 

One of the open problems of Hamiltonian perturbation theory is how to 
extend to infinite dimensional systems, at a finite specific energy (or tem- 
perature), the results known for systems with a finite number of degrees of 
freedom. Indeed there exist results both for infinite systems such as partial 
differential equation (see, for example, [HE]) or on infinite lattice systems 
(see OH]), but only for a finite total energy of the sistem, i.e., at zero 
temperature. 

In the present paper we provide perturbation estimates on the so called 
Klein Gordon lattice in the thermodynamic limit, at a finite temperature, 



*Universita di Milano, Dipartimento di Matematica, Via Saldini 50, 20133 Milano, 
Italy. E-mail: andrea.carati@unimi.it 



1 



by controlling, in place of the usual L°° norm, the L 2 norm relative to the 
Gibbs measure. If we denote by H the Hamiltonian of the system and by 
Ai the corresponding phase space, the Gibbs measure is defined by 



where Z(/3) = J M exp(—/3H(x))dx is the partition function, and /3 > the 
inverse temperature. 

We construct an adiabatic invariant whose time derivative has an L 2 
norm exponentially small in the perturbation parameters (see Theorem [T|). 
The construction of the adiabatic invariant is standard (see [5]), but the 
estimate of its time derivative in the L 2 norm involves some probabilistic 
techniques, which have been developed in the frame of statistical mechanics. 
In fact, since L 2 is not a Banach algebra, the usual scheme of perturbation 
estimates cannot be implemented. So the use of the algebra property is 
replaced here by a control of the decay of spatial correlations between the 
sites of the lattice, making use of techniques introduced by Dobrushin (see 
[6]). In particular, we are able to show that, for lattices in any dimension 
with finite range interaction (i.e., in which each particle interacts only with 
a finite numbers of neighbouring ones) the spatial correlations decay ex- 
ponentially fast with the distance. This requires also an estimate on the 
marginal probability densities induced by the measure [i on subsystems of 
finite size: this is done by adapting to lattices the techniques introduced by 
Bogolyubov et al. (see [7]) in interacting gas theory. 

The paper is organized as follows. The main result on the considered 
model, namely the construction of an adiabatic invariant in the thermody- 
namic limit, is stated in Section [2] (Theorem[T]), togheter with two corollaries 
concerning a control on the time evolution of the adiabatic invariant and a 
lower bound to its time autocorrelation. Then, in Section [3j we present the 
scheme of the proof of Theorem [IJ whereas the fundamental ingredients of 
the proof are separately given in the subsequent three sections. The first 
one (Section [4]) concerns perturbation techniques and deals with the for- 
mal construction of the adiabatic invariant. The other two sections have 
a probabilistic nature: the estimate of the marginal probability is given in 
Section [5]together with the estimate of the norm of the time derivative of the 
adiabatic invariant. In Section [HJ we state Theorem [2] in which the estimate 
of the spatial correlations is given, which enables us to give an estimate on 
the variance of the adiabatic invariant. The proof of Theorem [2] requires to 
apply a technique due to Dobrushin and Pechersky (see [8]), and is reported 
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in Appendix IB.3I In Section [7] we discuss how a lower bound on the time 
autocorrelation provides information on the relaxation time to equilibrium. 
The conclusions follow in Section [8l Most of the proofs of a more technical 
character are given in two appendices. 

2 Stabiliy estimate in the Klein Gordon lattice 

In the literature, as a prototype of several models, the so called Klein Gordon 
lattice is studied (see [9]-[l4]). Prom a physical point of view, it mimics a 
chain of particles, each free to move about a site of a lattice, subjected 
both to an on-site restoring nonlinear force and to a linear coupling with 
the nearest neighbours. It can also be seen as a discretization of the one- 
dimensional <£ 4 model, which plays a major role in field theory. 

The Hamiltonian of such a system, in suitably rescaled variables, can be 
written as H = Hq + H\ , in which 



where p = (pi, ■ ■ ■ ,Pn) and q = (<7i, ■ ■ ■ , <Zjv) are canonically conjugated 
variables in the phase space A4, and e is a positive parameter, while u is 

defined by ui d = \Jl + 2e. Since we don't want to face in this paper the 
problem of small divisors, which typically arises in perturbation theory, we 
confine ourselves to the case of small e, i.e, of small coupling between the 
sites. 

We aim at showing that, for small enough e and sufficiently large /3, 
there exists an adiabatic invariant for H (see Theorem Q] below). To come 
to a precise statement, we need some preliminaries. 

As usual, {X) will denote the mean value of a dynamical variable X with 
respect to the Gibbs measure fx relative to the given Hamiltonian H at a 
given (3, i.e., 





The L 2 (Ai,[i) norm of X is then \\X\\ = yj (X 2 ) and its variance a\ is 

defined according to a\ == {X 2 ) — (X) 2 . Finally, we also recall that the 
correlation coefficient of two dynamical variables X and Y is 




(3) 



PX,Y = 



axoy 
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and that X and Y are said to be uncorrelated if px,Y 
We can now state our main theorem, in which 
bracket, 



denotes Poisson 



Theorem 1 (Estimate on the adiabatic invariant) There exist positive 
constants e* , n, independent of N, such that if e < e* and j3 > e~ l , then 
there exists a polynomial function X uncorrelated with H such that 



[X,H] 
°x 



< 



exp 



1/4 



def 1 
~~ 1 



(4) 



Remark. We require X to be uncorrelated with H in order that our adi- 
abatic invariant be sufficiently different from the Hamiltonian, which is ob- 
viously a constant of motion. 

Before the proof, we point out immediately that this theorem has two rele- 
vant (and strictly related) consequences on the time evolution of the dynam- 
ical variable X. They will make clear in which sense i at the r.h.s. of (JH) 
can be seen as a stability time. The first consequence (Corollary [1]) concerns 
the probability P that the value of the variable X changes significantly from 
its original value. Indeed, it entails that the probability of such a change 
is practically negligible if t < i. The second consequence (Corollary [2]) is a 
lower bound on the time autocorrelation of X . We take here as definition 
of time autocorrelation of a dynamical variable the following one: 

Cx(t) d = px t ,x , 

where Xt(x) = X($*x), <&* is the flow generated by H and p the correlation 
coefficient defined by ([3]). We have chosen to rescale the usual definition, 
dividing it by because the variance of X is the natural scale of its 
autocorrelation, since Cx(0) = 1 and the inequality |Cx(£)| < 1 holds for 
any t. 

We report here both results, which follow from Theorem [T] and from the 
simple estimate \\X t - X\\ 2 < t 2 \\[X, H]\\ 2 . The latter can be found in the 
proof of Theorem 1 of paper |15| and is however reported here in Section [7] 
in order to make that section self contained. 

Corollary 1 In the hypotheses of Theorem^ for any A > one has 
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Corollary 2 In the hypotheses of Theorem^ one has 



*<«) > 



) 
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Remark. We observe that the notion of stability time for dynamical sys- 
tems is not unambiguously defined. In Section [7] we will provide a definition 
of "relaxation time" in terms of time autocorrelation of dynamical variables, 
which seems to us significant from a physical point of view. With such a 
definition, Theorem [1] turns out to mean that the "relaxation time" is ex- 
ponentially long in the perturbation parameters. 

3 Scheme of the proof of Theorem [L] 

First we use a variant of the classical construction scheme of approximate 
integrals of motion (see [16]) in order to perform the construction of the 
adiabatic invariant as a formal power series. Precisely, we use the scheme 
developed by Giorgilli and Galgani for a direct construction of integrals of 
motion (see [5] and Section|4]for the actual implementation). It is well known 
that the series thus obtained are, in general, divergent, so that the standard 
procedure consists in using as approximate integral of motion a truncation 
of the series. Denoting by Y n the series truncated at order 2n + 2, it turns 
out that it has the form 



where Pj are suitable polynomials. In order to make such a quantity un- 

def 

correlated with i?, it is convenient to consider X n = Y n — H instead of Y n 



In order to make the construction rigorous, one has to add rigorous 
estimates of the variance a\ of X n , and of the L 2 norm of [X n ,H]. The 
first step to get such estimates consists in controlling the structure of the 
polynomials Pj (which, in particular, contain only finite range couplings) 
and the size of their coefficients. This is done recursively, by a variant of 
the technique of the paper |17| . which is implemented in Section H] (see 
Lemma [2]). We emphasize that, at variance with the original paper, we 
obtain here estimates independent of the number of degrees of freedom. 

Then, due to the structure of the polynomials Pj, to get the needed I? 
estimates one has to compute the L 2 norm with respect to the Gibbs measure 



n 




(5) 



3=1 



itself. 
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of the monomials appearing in Pj . The key step for this computation consists 
in giving an upper bound independent of N to the marginal probabilities 
of the Gibbs measure. Such an estimate is obtained by adapting techniques 
developed by Bogolyubov and Ruelle (see |7j and [18]) and is reported in 
Lemma H] of Section [SJ One thus obtains the following bound 

X n \\ < ViV (vVf 1 (n!) 4 (/T 1 + e) n «? , (6) 

which is valid for a suitable constant K\ > 0, provided e is small enough and 
j3 large enough (see Lemma [3] of Section [5|) . 

We emphasize the presence of the factor yN and that K\ is independent 
of N. It will be shown that actually the l.h.s. of © is of order yN even 
if it is the square root of a sum of 0(N 2 ) terms. This is due to the fact 
that most of the terms have zero mean because the measure is even in p and 
furthermore the p's are independent variables. 

To get the Theorem, one also needs an estimate of ax n from below. This 
is obtained in two steps, which are based on the remark that ax„ > ®x x —ciii 

def 

where 1Z = X n — X± is a remainder. 

First we compute explicitly X\ and estimate from below ax 1 , obtaining 
a bound proportional to y/N. Then, we estimate from above an- Precisely, 
we use techniques introduced by Dobrushin in papers [U [8] to show that 
an behaves as \/N (see Lemma [8] of Section [6]) . We remark that this is the 
analogue of the law of large numbers. We recall that Dobrushin's techniques 
enable us to show that spatial correlations between variables pertaining to 
different lattice sites decrease exponentially with the distance between the 
sites, so that the monomials appearing in Pj are essentially independent, and 
the variance of P n is essentially the sum of the variances of each monomial. 
This leads to Lemma [7] of Section [61 which shows that, for small enough e 
and large enough (3, for n < k 2 1 ^ 4 (e + /3 -1 ) -1 / 4 there holds 

ax„>v^V(e + /3- 1 )/(8/3) , (7) 

where again K2 is a positive constant. 

Then one finds the optimal n, call it n, such that the ratio || [X^, H\ W/ax^ 
takes the minimal value. Notice that, as n belongs to a bounded domain, the 
minimum can be attained at the boundary. The optimization is immediately 
done, once the estimates are given both for the L 2 norm || [X n , of the 
time-derivative of the quasi integral of motion X n , and for its variance 
a\ . Then, the function X satisfying ([4]) of Theorem Q] is simply given by 

X d = X n - Hpx^HVXn/oH- The identity a\ = (1 - Px^u) a \^ together 



6 



with the upper bound to px n ,H given by Lemma enables us to extend all 
conclusions from X* to X. 



4 Construction of the adiabatic invariant 

Following |5], we look for the formal integral of motion by looking for a 
sequence of polynomials x = {Xs} s >\ such that 

[H, T x Hq] = at any order, (8) 

where T x is a linear operator, whose action on a polynomial function / is 
formally defined bj{3 

T x f = £ (T x f) s , with (T x /) = f / , (T x f) g = ^-IXj, (T x f) s _ 3 ] . 

s>0 j=l 

(9) 

Inserting the expansion of T x Hq and H in ([8]) and equating terms of equal 
order one gets the system 

Q Q = H , e s - L oXs = ^ s for s > , (10) 

where 



•'i r , . , i=i . (ii) 



= - E - [». (^o)._,] " E ( T x O^ for . > 2 , 
i=i z=i 

Lo = [-ffo, •] is the homological operator and (110p has to be read as an 
equation for the unknowns Xs, ©sj which have to belong, respectively, to 
the range and to the kernel of the operator Lq. By defining the projections 
ILv, fl^, respectively on the kernel Af and on the range 1Z of Lq, one thus 
determines recursively 

Xs = -L^ILr,*. , 9 S = U X ^s for s > 1 . (12) 

The approximate integral of motion is then obtained by truncating the se- 
quence T x Hq at a suitable order. 



1 Notice that in paper [5] the \s were required to be homogeneous polynomials of degree 
s + 2. However, there is no problem in considering the present more general case. 
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We have to estimate the action of the operator T x on the class of func- 
tions f(p,q) we are interested in, in a norm which is well suited for our 
problem. Such a norm is defined as follows. Let Hi' 1 denote the class of 
monomials^] p k q l of degree s, i.e., with + \l\ = s, which furthermore de- 
pend on sites that are at most r lattice steps away from i, namely such that 
kj = lj = if \i — j\ > r. We denote by V s ,r the set of all homogeneous 
polynomials of degree s that can be decomposed as 



N 



with fij G Us' 1 , where 



associate a normjf) defined by 



is the cardinality of T-Q 1 . To / G V s r we 



def . 
= mm < 




max > leal > , (14) 
ie{l,...,N} p[ 3 ( 



where the minimum is taken over all possible decompositions of /. 

Now, we can estimate the action of T x on any function / G V s ,r according 
to the following Lemma, which is proved in Appendix 1X1 

Lemma 1 LetT x be the operator defined by relative to the sequence x = 
{Xs}s>o which solves the system of equations (E2H22P for the Hamiltonian 
HI). Then, for any f{p,q) G V2 S +2,r, one has (T x f) n = YZ=Q fn > wher e 

fn +l ^ G V2s+2l+2,r+n-l and 



■J n 



(s+l) 



" 1 1 n - r\ s\ l\{n-l)f 

(15) 



Lemma [2] below will give bounds to the adiabatic invariant obtained by 
truncating at a finite order the formal power series which defines T x Hq. In 
particular, the adiabatic invariant will simply be Y n = Yl™=o (Tx-^o). > so 
that the polynomials Pj appearing at the r.h.s. of ([5]) of Theorem [T] are 

Pi * (T X H ) 3 , (16) 



2 We adopt here the multi-index notation: k = ki, . . . , fev and I — h, . . . , In are vectors 
of integers, with \k\ = \ki\ + . . . + \k 2 \. So, p k q l = p^ 1 ■ . . . • p% N q[ x ■ ... • gjf . 
3 One can check that this is indeed a norm. 



s 



while the quantity we will focus on will be 

n 

X n = Y n -H = -Q l + J2 p j 

3=2 

The time derivative of X n is then given by 

X n = [X n ,H] = [P n ,Hi] , 



(17) 



(18) 



which is a polynomial of order 2n + 4. In order to obtain the estimates 
of the L 2 -norm, eventually, it is of interest to take into account the parity 
properties of the operator T x , with respect to the canonical coordinate p. 
So we define as V + the space of polynomials of even order in p, and V the 
space of those of odd order in p. 
Finally, we can state 

Lemma 2 For the adiabatic invariant constructed through T x Hq (see Mb}) ) 
one can write 

n . 

Pn = Y. 777^7^4° , (19) 



^ l\(n-l)\ 



where p}p £ V + Pi V 2 i+2,n-i and 



p(0 

± n 

Furthermore, one has 



< V n , with V n = T- i n\ 



n+l 



(20) 



with x [ n ] ern V21+2 



< C n , with C n = 48 • 2 12n n! ((n + l)!) 2 . 



(21) 



Proof. The proof of the upper bounds is mainly based on the application 
of Lemma Q] to the function Hq G V201 together with the simple bound 
|| JEZq ][■ = uj < 2, which holds for small enough e. This proves equations 
(fT9|) . (|20|) . Then, we use the fact that [X n ,H] = [P n ,H{\ and the upper 
bound to the norm of the Poisson brackets of two variables provided by 
Lemma [TU1 of Appendix |A| This gives equation (|2ip . 

The parity properties are obtained by observing that [V^, V^] C V + 
and |'P ± ,'P =F ] C V~, as well as U A f(V + ) C V + and ILyCP") C V~ and that 
the similar inclusions regarding n^hold, and then working recursively. 
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Q.E.D. 



5 Marginal probability estimates 

The aim of this section is to prove the bound on the norm of X n given by 
the following 

Lemma 3 There exist constants j3 > 0, e > 0, k\ > such that, for any 
f3 > /3 and for any e < e, for X n defined by hlSfy of Section [7] one has 



X n 



< VN(V2p] \n\) A {p- 1 + e) n ^ . (22) 



The key tool of the proof is an estimate of the probability that the co- 
ordinates of a finite number s of sites are near some fixed values. Such 
an estimate is given in the following Subection 15.11 whereas the proof of 
Lemma [3] is given in Subection 15.21 

5.1 Estimates on the marginal probability 

Everything is trivial for the p coordinates, for which the measure can be 
decomposed as a product: from a probabilistic point of view, this means 
that every pj is independent of the q and of any pi, for i ^ j. We focus, 
instead, on the q coordinates, which are independent of the p, but depend 
on each other. Then, we must study the relevant part of the density, which 
is given by 

D N (q 1} . . . ,q N ) d = — exp[-/3U N (q 1 ,...,q N )] , (23) 

where Zn is the "spatial" partition function 

def f +ao r+oo 
Z N = dqi... dq N exp [-(3U N (qi, q N )\ (24) 

J— oo J — oo 

and Un the part of Hamiltonian ([2]) which depends on q, namely, the po- 
tential 

N / 2 a \ 2V-1 

def / qf , qf \ , UU+i 



iv , 2 4 \ iv-i 

i/*ta.... W )*E(«£+&)+«E „ 

i=l v ' i=l 



The main point is then to estimate the marginal probability Fg^P (q^ 
that we are going to define. Given a set of indices i\ < 12 < ■ ■ ■ < i s we 
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say that they form a connected block if ij + i = ij + 1, i.e., if they label a 
"connected" chain. We say that a sequence of indices i\ < %2 < ■ ■ ■ < i s 
form y blocks if the set {ij}j =1 can be decomposed into y connected blocks, 
which furthermore are not connected to each other. Given a set of indices 
ii < i 2 < ■ ■ ■ < i s we define 



(N) def [ + °° f + °° 

Fs, t (%i, ■ • • , 1i s ) = / dq is+1 ... dq iN D N (q u ...,q N ) 

J -co J— CO 



(25) 



where y is the number of blocks in the set {ij} S j=\- We remark here that 
such a quantity depends on the number of particles, N, but we will find for 
it an upper bound independent of N. In fact, the estimate will depend only 
on s and y, but not on the precise choice of the sites. 
Define the two functions 



def 



n s Mhi- ■ ■ ,Qi a ) = ex P 



C \ . ^ c (Qi k -QiiY 



\k=l \ / k,l=l 



2uj 



V k=i 



n s , f (qh,---,qi 3 ) = exp 



where <5,- is the Kronecker delta. 




CO 



(26) 



(27) 



Remark. Notice that n SiJ is the configurational part of the Gibbs measure 
of the system with variables % , . . . , qi s and free boundary conditions, apart 
from the absence of the normalization factor (i.e., the partition function), 
whereas h SJ is the analogous quantity for the same system, but with fixed 
boundary conditions. Thus, they differ only because of the different depen- 
dence on the coordinates at the sites lying on the boundary of the blocks, 
the number of which, 7, satisfies y < 7 < 2y. If we denote by mi, . . . , m 7 
the indices of these sites, we can write the identity 



n s,z{qh • 



JJexp 



(3e 



-a 



< 



j=i 



(28) 



where the factor a mj is equal to 1 or 1/2 according to whether the site rrij 
is isolated (i.e., the block is composed of only that site) or not. 
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Then the following lemma, which is the main result of the present sub- 
section, holds 



Lemma 4 There exist constants /3 > 0, e>0 ; K > and a sequence 
> such that, for any j3 > (3 and for any £ < e, one has the inequalities 




(30) 

The proof of such a lemma is based on the techniques of paper [7] , which 
apply quite simply to the case of periodic boundary conditions (see Lemma[5] 
below), on account of the translational invariance. Thus, it is also useful to 
introduce the density Dn relative to the periodic system, defined by 

def 1 

D N (qi, ...,q N ) = — — exp [~pU N (qi, ...,q N )+ f3eqiq N ] . (31) 

In this definition there appears the partition function for the periodic case 

dcf / ,+ °° r+oo 
Q N = dqi... dq N exp [-/3U N (qi, ...,q N )+ fieqiq N } . (32) 

J —oo J —oo 

For the periodic system it is simple to estimate two relevant quantities. 
The former is the ratio between the partition function for N — 1 particles 
and that for N particles, i.e., the ratio Qn-i/Qn- The relation between 
Qn and Zj^ is then obtained as a particular case of Lemma [6l which will be 
stated later on. The latter is the probability, evaluated with respect to the 
density for iV particles, that the coordinates of r particles have an absolute 
value smaller than Oy2w//3, for a given 0. In other terms, we need an 
estimate of the following quantity 

A f2u /2u7\ dcf f +co f + °° 

Pistol < ©V J A • • • A < V TJ = Lj qi ■ ■ -J_J qNl M<eV^ 

x ••• x Wev^^ 1 '---'^ ' 

(33) 

in which 1a is the indicator function of the set A. We can now give the 
mentioned estimates by the following lemma, whose proof is deferred to 
Appendix IB. 1[ 
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Lemma 5 There exist constants (3q > 0, Eq > 0, Kq > 2 such that, for any 
(3 > (3q and e < Eq, one has 



Qn-i „ „ P 



q^H^L- (34) 



Furthermore, if > 2^/r \og(ArKo), one has 

Pjv(|gi| <@^A...A\q r \ <Q ]fj) ^ \ ■ ( 35 ) 
This result enables us to give the proof of Lemma HI 

Proof of Lemma [4] First write 

D N (q 1 ....,q N ) = n N _ S j> (q is+1 ,...,q iN )) n s ^ (q h , . . . , q ia ) I (q u . . . , q N ) , 

for a suitable p', with p — 1 < p' < f + 1 (the lower and the upper bound are 
attained, respectively, if both 1 and N are contained in ii,...,i s or none 
of them), where / contains the terms of interaction between the "internal" 
and the external part of the system. Remarking that I < 1, one gets 

. . ^ / r+oo r+oo 



N 



x n N - S ji(q ia+1 ,. . .,q iN fj n S j(q h ,. ..,q is ) .(36) 



We now have to estimate the integral appearing in (|36p . More in general, 
in the course of the proof we need to estimate integrals of a similar type. 
This will be done in Lemma El which will be given in a while. Introduce the 
quantities: 

-oo r+oo 



def r+°° r+oo 

Q l M = inf / dq 1 ... dq M n M Aqi, ■ <?m) 



def f 

Qm = SU P / d <?i 

B(M,f) J-oo 



oo r+oo 



(37) 



where B(M,jc) denotes the collection of all possible partitions of M indices 
in p blocks. It is also convenient to consider the quantities defined in a 
similar way, by integrating fiM,i in place of um,i, namely 

'•+oo r+oo 



def r+°° r+oo 

Z\, = inf / Mi-- - I dq M n M ,p(qi, qhi) , 

def r+°° r+oo 

Z l M = sup / dqi... dq M n M ,i{qi, ■ ■ ■ ,qu) ■ 

B(M,t) J-oo J-oo 
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It is easily shown that for p = 1 one has Z M = Z M = Zjy. In order to link 
them to Qn and to each other, we use the following lemma, the proof of 
which is deferred to Appendix IB.2I 

Lemma 6 Let /3q > 0, Eq > and Kq > 2 be constants such that Lemma [5| 
holds. Then, for any f3 > /3q and any e < Eq, the inequalities 

^>1, j± > \ (8tK )- 32 ^ 2 (39) 
hold. Furthermore, the chain of inequalities 

^<®k<2Kl exp (4 P£o k( ? , K q )) , (40) 
holds, where R(p, Kq) is the solution of the equation 

K^n h R) = \, (41) 
T(s,x) being the upper regularized Gamma function 

j^wC^ 1 '"*- (42) 

The previous lemma enables one to see that Z^ 1 < 2(8K ) 32e °/Q N , while 
the integral appearing in (j36|) is estimated by 

" ^N-s 



Q* N _ s < 2K* +l exp [4(y + l)e «(? + 1, Kq)} Q 
Thus, due to relation ([34"j) of Lemma [5l one easily sees that 

so that ([29]) is proved, taking 

£ ? > 2 9fe »+ 2 ^ +1+32eo exp [4(y + l)e «(? + 1, K )] . (44) 
We come now to the proof of (|30p . To this end, we write 



Qn-s jj 
Zn 

x G(q mi ,. ..,q mj ,q h ,.. . ,q, ,) 



• • • ,<?jv) = — D N - S (q is+1 ,. . . ,qi N )n S j(q h , . ..,q is ) 
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where the sites Zj are the ones which are contiguous to the blocks, but not 
contained in them, taken by keeping the relative order. Furthermore, due 
to the periodicity there appear factors depending on q\ and qw, if the sites 
1 and ./V are not contained in i\, . . . , i s . In this case, we put = q± 
and qi , = qjy. We denote by 7', with 7' < 2r + 2, the number of such 
indices. The explicit expression of the function G is complicated. Plainly, it 
represents the product of the factors exp(— Y2j P^QmjQk/ 1 ^) among all sites 
Zj contiguous to nij, and rust the factor exp(/3eg/ i g , ; i+1 /a;), when Zj and Zj+i 
belong to different blocksu In any case, a lower bound to G in the region 

A = <2 ?v / log(4K )v / 2^ A ... A | fty |<2y v / log(4K'c 

is given by 



G > exp I -4^A^(4^)^ \ qm .\j (4K )- 8 ^ 2 +f) . 

So, we can write 

FljPfa,. ..,«.)> ".(«x,. • • ,0(4tfo)- &(?2+r) ^PiV- a (-4) 

x exp -4erW ^0og(4K o )^ \q. 

where the (positive) contribution of the integral over A c was neglected. 

The term with Pn_ s (A) in (j4"5j) is bounded from below by relation (f35l) 
of Lemma [5] As for the fraction, by Lemma [6] we obtain 

1 

2Koexp(8FoR(l,Ko)) 

Now, operating as in the deduction of formula (03]), it is sufficient to observe 
that Q] v _ 1 ^ v //3/(2vra;)Q] v to obtain Q 1 N _ S > ((3 /2(mj)) s/2 Q l N . Then, 
choosing e, /3 such that Kq < e 4 /4 and observing that Q} N > Zjy, one gets 
(130]) with 

£ r > (4^ ) 8eo{r2+r)+1 exp(fe «(l, #o) • (46) 
Finally, <t r can be chosen as the maximum of the r.h.s. of (|44p and of (|46p . 
This concludes the proof. 

Q.E.D. 



Qn-s > nT ^ ; , _ ^ r ^ ^ Qn-s 



4 Notice that the index i lies in {1, . . . ,7'}. But if in the expression exp(/3eqi i qi i+1 /u 
there appears qi , then one has to intend simply £ 7 /+i as h. 
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5.2 Estimate of \\X„ 



We apply directly inequality (f2Uj) to get the proof of Lemma EJ using the fact 
that such a quantity is a sum of polynomials depending at most on 2n + 3 
sites, as can be seen by Lemma [2] of Section HI 

Proof of Lemma[3]The key ingredient of the proof is, as stated in SectionHl 
that the polynomials P n are even in the p coordinates, so that the X^s are 
odd in the p. On account of that, X^ is a sum in which the terms coming 
from the product of two monomials depending on separated groups of sites 
contain at least one pi to an odd power. Since the measure is even with 
respect to any p, these terms have a vanishing integral. 

We formalize this way of reasoning by decomposing X n as X n 
where the /j's are polynomials depending at most on the sites between i — 
n — 1 and i + n + 1. Then, the L 2 -norm of X n is expressed according to 



Si=l fii 



X r . 



~^2dj=i{fifj) • I n this sum, all the terms with > 2n+2 vanish, 



while the other ones are estimated in terms of 
On account of Lemma [21 we can write 



X r . 



in the following way. 



n+l 



iT/n+l— l,i I 
\ n 2l+2 I 



_ ^ (n + l)! ^ n+1 _; (o 
h ~ 2^ l\( n + l-l)\ ^ ls ' lhs ' 

1=0 y ' s=l 



in which f- l J is a monomial in H^+T^ anc ^ ^ e decomposition in these 
monomials is performed in such a way that sup^ |cj Sj /| < C n . Then, we 
sum on j and obtain that 



X r . 



N n+l 



<(4„ + 5 ) C ;|EE il ,'!" + ^n, ^ +2 -' W. .>>■ 



l\(2n + 2-l)\ 
=i i=o v ; 



geH 



2n+2-l,i 
21+4 



where we used the fact that the only nonvanishing contributions to the 

S- r ) ^ ijn+l-l+rj i Ar) ^ nl n+l-r,k 
^2«-2r+2 dnu J km 



integral come from the product of / - s G ^2l-2r+2 ano - J km. e ^2r+2 

def Al-r) Jr) qJ 2n+2-l,i 



for \j - k\ < 2n + 2 - I, so that g = 
between j and /c. 



fkm ^ ^ 



2/+4 



, for a suitable z 



Then, we make use of (I29p together with the estimate (126|) for n s ^ to 
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bound the mean value of any function in T~L^i^ In fact, one has 



SUp (J^^JI I" : 

,,2n+2-l,i V Z7TW J—oo 



9^21+4. 



„. ,^'V (2/ + 3)! 

1 1 /3 J 2 l + 2 



x 2l+4 exp ( ——x 2 ) dx 

Z+2 



So, the inequality 

2 / 1 \ 2ra + 2 

< CiET 4n+5 (4n + 5)(2n + 4)! (2w) 2n+4 /T 2 U + 

holds. Thus, choosing a suitable «i > and using the value of C n given by 
(|2ip of Lemma [21 inequality (122j) is satisfied. 

Q.E.D. 



6 Estimate of the variance of the adiabatic invari- 
ant 



In the present Section we prove the following Lemma which was used in 

4 



the proof of Theorem [TJ The lemma concerns estimates on the variance a 2 



and on the correlation px n ,H of the adiabatic invariant and reads 

Lemma 7 There exist positive constants e > 0, K2 > 0, K3 > 1, such that, 
for any e < e, for any (3 > and for n < k 2 1 ^ 4 (e + (3~ 1 )~ 1 / 4 , with X n 
defined by |i7| j ; the following inequalities hold: 

ax n >^N £ -^f- (47) 

and 

l 1 £ 2 V"" 2 

The proof of this lemma requires the study the spatial correlations be- 
tween quantities depending on two separate blocks. The study of these 
properties has to be performed within the general frame of Gibbsian fields 
and conditional probabilities. In the present Section we provide the neces- 
sary notions and give a proposition of a general character concerning the 
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decay of spatial correlations for lattices with finite range of interaction, i.e, 
Theorem [2] of Subsection I6.2[ from which it will be possible to finally come 
to the proof of Lemma which will be given in Subsection 16, 4i 

Our treatment of conditional probabilities is inspired in particular by the 
work of Dobrushin (see [B]). More precisely we will make reference to paper 
[6] for the main ideas, and to the subsequent beautiful but underestimated 
subsequent paper [S], by Dobrushin and Pechersky, for a more direct relation 
to our problem. As a matter of fact, most of the ideas of this section are 
already contained in works [6] and [8], but the explicit result on the spatial 
correlations given here required some additional work. We recall that, since 
Gibbsian fields and the related techniques were introduced in order to deal 
with infinite lattices, our result holds even if the number of sites tends to 
infinity. 

The present section is structured as follows: in Subsection 16.11 the link 
between spatial correlations and conditional probabilities is shown, and in 
Subsection [62] we state Theorem^ whose proof is deferred to Appendix lB.3l 
Such a result is used in order to obtain an upper bound to ap n , stated in 
Lemma [8] of Subection 16.31 whence the proof of Lemma [7J easily follows, as 
shown in Subection 16.41 

6.1 Link between spatial correlations and conditional prob- 
ability 

In order to prove Lemma [7] we have to estimate quantities such as (fg) — 
(f)(g), relative to the Gibbs measure fj,, where / is a function which depends 
on sites belonging to a set V, while g depends only on sites in V, with 
V PI V = 0. Our aim is to show that such correlations decrease as the 
distance between V and V increases, where the distance d(V, V) is defined 

for example as d(V, V) = f inf ig y \i — j\. 

We start showing the relation between the spatial correlations and the 
conditional probability in a setting more general than ours. We consider as 
given a measure fj. on with TcZ", which induces on the measurable 
set i C R 1 ^ 1 the probability 

p v( A ) d = / dM^J-Axnv-O) . 

Jr\t\ M 

where 1a is the indicator function of the set A. One can express the quantity 



18 



we are interested in as 

(fg)-(f)(g) = [ J(x)Pv(dx)([ s(y)iV(dy|dx)-/ g(y)P v (dy) 

(49) 

where Py(B\A) represents the conditional probability of the measurable set 
once A is given. So, in order to estimate the correlation between 
two functions, it is sufficient to estimate the difference enclosed in brackets 



at the r.h.s. of (|49p . Now, we notice that for any pair of probabilities P and 
P on , one has 

/ g(x)P(dx)- [ s(y)P(dy) < f (\g(x)\ + \g(y)\) 

JrW\ JrW\ Jr\v\ x r\v\ (50) 

x l x ^ y <5(dx,dy) , 

in which Q is any probability on M) v \ x Rl^l such that P and P are its 
marginal probabilities. In other terms Q is a joint probability of P and P, 
i.e., for any measurable B C W v \ one has 

P{B) = Q(Bx Rl y l) and P{B) = Q (m) v \ x b) . (51) 



Remark here that Q is not unique: indeed, such a probability provides also 
a way to define a distance between two probabilities defined on the same set 
V of indices, by 

D(P, P) d ^ f inf / l^y Q(dx, dy) . (52) 

We stress that the infimum is attained, i.e., there exists a probability mea- 
sure Q(dx, dy) such that D(P,P) = J l x ^ y <5(dx, dy) (see Lemma 1 of pa- 
per [8]). For the following, we suppose that it is possible to find a compact 
function /ijl with domain in R, such that 

\g(x)\ <Y,h(xi) , (53) 

i£V 

as is the case for the monomials we are dealing with. The bounds will then 
be given in terms of h. Now, observing that l x ^y < SieV ^-Xi^yn we can 
rewrite (1501) as 



f </(x)P(dx)- I g{y)P(dy) < £ I (h( Xj ) + %,)) 
■/aivi ymvi ij^t/^ 2 (54) 



X 



W yi Q(dx,dy) . 



5 See later for the definition of a compact function, according to the convention of paper 

0. 
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This way, we can make a direct connectino with paper [8], in which the 
problem of estimating the r.h.s. of the above expression is dealt with. We 
summarize here the results and the methods we need. 

6.2 Main argument and theorem on correlations 

In the quoted work [8], the framework is more general than ours, because 
it deals with the problem of defining a "probability" for the configuration 
of an actually infinite z^-dimensional lattice of particle, in terms of the set 
of conditional probabilities on each site, which is called the specification 
T. Now, our case is in principle different, because our lattice is finite, and 
the probability is defined through the Gibbs measure. In particular, the 
specification too is assigned by such a measure. 

However, as proved in [8], under suitable assumptions assigning the spec- 
ification uniquely determines the probability, i.e. the Gibbsian field, which, 
in our case, turns out to be precisely that of Gibbs. So, in our case, it is 
equivalent to speak in terms of specification or in terms of measure. Indeed, 
in this subsection we will speak in terms of specifications, and in the follow- 
ing one we will show that the specification determined by Gibbs measure ([I]) 
with the Hamiltonian ([2]) satisfies the assumptions of [8] (i.e. Conditions [T] 
and [2] below). 

We notice that the r.h.s. of (|54p can be bounded from above if one 
estimates the quantity 

V^Mtf)]} . ( 55 ) 

where £ and £ 2 are two Gibbsian fields which assign, respectively, the prob- 
abilities P and P appearing in (I54p and the expectations are obtained by 
integrating over a joint probability Q of the two fields. Indeed, in [8] an up- 
per bound just to A(j,i) is given, by requiring that two suitable conditions 
are satisfied. 

So, by adopting the same techniques of [8] we bound from above the 
r.h.s. (I54D and, thus, the l.h.s. of (I49D . Such a bound is given in Theorem [2] 
below. In order to state it, we recall the main notations of [8]. 

First, we consider a lattice of sites contained in T C Z", and a finite- 
range specification with a radius of interaction r (this means that the condi- 
tional probability at site i does not depend on the conditioning at sites j for 
|i— j| > r). Then, for a vector x G RI t I we denote by Pi tX (dx) the probability 
distribution conditioned to x everywhere but at site i. The specification T 



A(j> i) = max { E 



,E 
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is defined by 




Furthermore, we will say that a continuous positive function h on a metric 
space X is compact if, for any k > 0, the set {x G X : h{x) < k} is compact. 

For a fixed integer r, let <9 r V = f { j G T : j V, minify |j — k| < r} be 
the boundary of thickness r of a set V C T. We call a the number of indices 
such that |i| < r, i 7^ 0, where r is the range of interaction. 

If Zq is a maximal subgroup of Z 17 satisfying the condition |j — k| > r, 
for j,k G Zq, we denote by b the number of elements in the factor group 



The conditions of paper [8] (which are hypotheses on the specification 
r, once the compact function h is given) are the following 

Condition 1 (Compactness) Let h be a compact function on R and let 
C > and cj > 0, for |j| < r, j 7^ 0, 6e some constants. We suppose that 



Condition 2 (Contractivity) Zei if > and k- } = kj{K) > 0, for |j| < r, 

j 7^ 0, be constants and h be a compact function. We suppose that 



TL v \Zq. 




|j|<r,j#) 



2. for any i G T and any x G R' T ' one has 



h{x)P i>x (dx)<C + Cj_i^(xj) 
je9 r {i} 



def 

. a = 



|j|<r,j^0 



/or any i G T and any pair 0/ configurations 



x\x 2 G K |T| sitc/i f/iai 



je9r{i} 



max 




one /ias i/ie inequality 



jeO r {i} 



where D(-,-) is the distance defined by {52^). 
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The set of specifications (i.e., the sets of conditional probabilities on 
every site) which satisfy Condition [1] for the constants C, 5 and the com- 
pact function h will be denoted by @(h,C,5). We will instead denote by 
A(h, K, a) the set of specifications which satisfy Condition [2] for the con- 
stants K, a and the compact function h. 

If the specification satisfies Conditions [T] and El Dobrushin and Pech- 
ersky show that the specification uniquely determines the probability (see 
Theorem 1 of [8]). In particular, also the marginal probability Py on V, and 
the probability Py(-|dx) conditioned to a vector x in V are determined by 
r, and the maximum value taken by A(j,i) of (|55p is bounded from above. 

As a matter of fact, in [8] the authors do not investigate the explicit 
dependence of A(j, i) on |i— j|, which is what we need, but such a dependence 
can be obtained by a slight modification of their way of reasoning, which 
leads to the proof of Theorem [2] that appears in Appendix IB. 31 We need 
also to introduce an upper bound to the mean value of h with respect to the 
probability Py on V and to the probability JV(-,dx), i.e., 

Wx = sup max! f %)fl(dj/|dx), f h(y)I\(dy),c\ , (56) 
i&V Ur Jr J 

where x E V, while C is the constant entering Condition [TJ So we can state 



Theorem 2 (Decay of Correlations) Let f be a measurable function from 
Rl^l to R, depending on the sites lying in the set V. Let g be a measurable 
function from R^ to R, depending on the sites contained in the set V , with 
VC\V = 0, and h be a compact function such that inequality A53\) is satisfied. 

Let r E @(h,C,S). Then, there exists a constant Kq, depending on 
h,C,a,b only, such that, if T S A(h, KqC,ol) n @(h,C,5) with any a, one 
can find constants D,c > 0, for which one has 



\(fg)-(f)(g)\<D\V\ 



I V| 



f( x ){h)Mdx) 



exp ( -cd(V, V) ) . (57) 



The constant D depends on a,b,a,8 only, while one has 



def 1 . 

c = -7- log 
or 



1 



max{a, 5a b } + 1 



(58) 



6.3 Estimate of the variance of P n 

The previous way of proceeding can be fitted to our case by choosing as 
specification that given by the Gibbs measure relative to H. As / and g, 
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we choose polynomials in p and q, depending on two disjoint sets of sites. 
In fact, on account of Lemma [21 we know that the P n are constituted by a 
sum of such terms. This way we can study o~p n and state that it is bounded 
from above as in the following 

Lemma 8 There exist constants f3 > 0, e > 0, k' > such that, for any 
j3 > j3 and any e < e, one has, for n < 1/e, 

a Pn <VNn\(k') n V n (V2py 1 (e + rT , (59) 
where the polynomials P n are defined by 1116}) and T> n are given in Lemma\^ 

Proof. Lemma [2] provides the necessary estimates for the coefficients which 
appear in the sum defining P n (see (fTM2"0j) l: we can write P n = X^=i/«j 
where fi are polynomials depending at most on the sites between i — n and 
i + n. The variance can be expressed as ap n = Ylij=i (ififj) ~ {fi)(fj}) ■ 

We then consider the set Si = f {(i,j) ■ \j — i\ < 2n} and ^2 = Sf. The 
proof goes on by finding an upper bound separately for the contributions 
coming from these two sets: for the latter, we use the methods developed in 
the present section, while the terms of the former group are estimated in a 
way similar to that of Lemma El 

We start from <Si. Firstly, we observe that, in general, one has 

\ififj) - (fi)(fj)\ < o U a h < max{</2>, (//)} ; 

so that it suffices to evaluate supj(/,f) in a way similar to Lemma [3] and to 
sum over the j with \j — i\ < 2n to see that the contribution due to the 
terms in this set is smaller than 

E KM') -</*></;>! <^K in+ \An+l)(2n+l)\{2u) 2n r 2 L + -X l NVl. 

We come now to S2- We will show below that the specification coming 
from the Gibbs measure satisfies the hypotheses of Theorem [21 and we make 
use of it in estimating the terms in the set 52 in the following way. We 
separate the terms of a definite degree by writing 

\%l n — hi I 

n . \ n 2l+2 I 

f. = V— e n ~ l V c ,f il) 

Z=0 V ; 8=1 
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in which f>® is a monomial in H21+2 anc ^ sn Pi,l zCs \ c is,l\ — ^n- We fix an 
index i and use Theorem [2] with f = fi and g = /jp, for every j / i; then, 
we sum over Z, s, j and i, subsequently. For each / we choose the compact 
function hi(x) as |x| 2 ' +2 , which satisfies (|53|) for any /j^. We will show that, 
for /? large enough and e small enough, one has 



(60) 



for a suitable constant k. Then the sum on s brings in a factor D n . As 
regards the integration in V, we observe that we can write 



P 2 1 8£ V2 W lj| -4 J 



1 



provided e is small enough. Thus, there exists k > such that 

2n 



\m)-{fi){m<~k n (n\fvl 



1 



e + 



exp [— c(|i 



2n- 1] 



where the constant c is defined by ([58 j) . and is, in the present case, equal 
to log(4/3)/2. Since the sum over j of such terms converges as N — > 00, 
the proof will be concluded if we show that (|60p and the hypotheses of 
Theorem [2] are satisfied. 

We proceed as in the proof of Theorem 2 of paper [8], starting from the 
explicit form of the conditional probability distribution given by the Gibbs 
measure: one has 



1 



exp 



-P 



2 4 

wr x e e 

—Z- + -7—? + ~ X • X i~l + ~ X ■ X i+l 

2 4or 



in which there appears the conditional partition function 

ry def 



exp 



2 4 

ux x e e 

+ -7—0 + ~ x • x i-l + ~ x • x i+l 

2 4w^ 



a; 



dx . 



As regards Condition [TJ we consider e < 2 5 fixed and define 

y(x) d = max{|xi_i| , \x i+ i\} and y(x) = f min|l/y(x), y^} • 

So, it is easily proved that, for j3 > 1, inequality Z x > cy(x)//3, holds for 
some c > indipendent of x, /3 and e. To prove it, it is sufficient to observe 
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that the integrand of Z x is bounded away from zero if \x\ < y(x)//3, and 
then to integrate over such an interval. We now show item 2 of Condition [1] 
for h(x) = \x\ 2l+2 . We note that 



/ hi(x)P i;X (x)dx < — x 2/+2 exp (-(3— )dx 

J\x\>y( X )/4 Jl/(4y( X )) V 4 / 

cy(x) 7v?/(4|/(x)) 

which, in turn, is smaller than (J + for a suitable constant -B, 

independent of e, f3 and Z, since the integral decreases as an exponential 
function of y/j3/y(x). Here, we have chosen (I + 1)! so that the previous 
relation is satisfied independently of I. Furthermore, one has 

/ h l (x)P t , x (x)dx<2- il - i (y(x)) 2l+2 
J\x\<y(x)/4 

and (hi(xi-i) + hi(xi+i))/hi(y(x)) > 1, for any x: this implies that 
J hi(x)P i>x (dx) <(l + 1)! I — J + —hi + — hj (a?i+i) . 



So, item 2 of Condition Q] holds with C = (I + 1)! {B/f3) l+l and d = c 
1/16. Since we have a = 2, b = 2, r = 1, item 1 of Condition [T] holds with 
<5a b = 1/2. 

Condition [2] is proved by computing two limits, first letting /3 tend to in- 
finity and then letting e — > 0. In fact, let x^ m \ for m = 1,2, be two different 
configuration such that < K / e^/ Bl / (e(3) and |ic^| < K /e^ Bl/(ef3). 

Then it is easily checked that 



-i 



lim 

/3— >oo , 



/ ^ x(1) (x)-P ix(2) (x) dx= / dze- z2 f(e,z, Z W)-f(e,z, Z ^ 

JR JR 



where 



,(m)\ def / _£ / (m) (m)\ / (m) (m) N 2 



CO 



f(e, z, *W ) ^ exp ( --z [ zH + z%> ) + — 2 [ + z, 



-i+i 



and z = x\/~j3, 2^ = x^y/jB. Now, by the dominated convergence 



j 



theorem, one has that the limit for e — )■ of f(e,z,z^) is equal to 1. 
Here, use is made of the fact that | can be bounded from above by 
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eK I e\J Bl I e < K je\J Be/e, because I < n, and n is smaller than 1/e, by 
hypothesis. So, for (3 sufficiently large and e small enough one has 



/ \ p i,^(x) ~ Pi,^(x)\ dx < - . 

JR 4 



We have chosen a bound to x m of this particular form, because the constant 
K in Condition [2] turns out to be smaller than K 2l+2 , so that it is inde- 
pendent of /?. So Condition [2] holds with K = K 2l+2 , k-\ = h = 1/2 and 
a = 1/2. Thus, Theorem [2] holds. 

There still remains to estimate {hi}*, in our case. By looking at its defi- 
nition (|56p . we notice that we have to estimate the integrals J h(y)P{(dy\dx) 
and J h(y)P\{dy). Now, on account of LemmaHJand relation (|28p . the dis- 
tribution functions of Pi(dy|dx) and P\{dy) can be bounded by 



e -Py 2 /{2u) 



Then, we use the bound to C previously found, together with the fact that 



2vrw J R y ' \P J 2 Z +! ' 

and we get ([60]) . This concludes the proof. 

Q.E.D. 



6.4 Estimate of the variance of X n 

Lemma [8] of Section [6.31 enables us to bound from below the variance of X n 
defined by (fTT|) and to estimate the correlation coefficient px n ,H, according 
to Lemma [3 which we prove here. 

Proof of Lemma [7] We start by recalling that, on account of (|17p . one has 
X n = —©i + $^? =2 -fjj with ©i defined by equations (fT2"HTTj) of Section HI 
It is easily seen that Bi = F + G + Hi, in which 

e N ~ 1 3 N 

i=l i=l 
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and IZi is the remainder. Then, we study the properties of F and G, for 
which the mean value, the variance and the correlation with H can be com- 
puted almost exactly, and we extend such properties to Oi, and to the whole 
X n , by observing that, in some sense, 0i is the term of first order in e + 

As regards formula (|47|) . we notice that, since F is odd in the momenta, 
while G, IZi and the measure are even, then F is uncorrelated both with G 
and with TZ±. Furthermore, one can observe that 

q N 

(GRt) - (G) (Rt) = -jjj-j «pj) - ( P i)(p*)) , 

and use the estimates of Lemma [5] to bound from above (qf), in order to 
prove that cr^ > a 2 F + a% + 2C GjKl > N(e 2 + /3~ 2 )/(8/3 2 ), where the second 
inequality holds for e and /3 _1 small enough. On the other hand, making 
use of ([59]) of Lemma [8] together with the estimate for V n given by ([20]) of 
Lemma [21 one has 



for a suitable constant K2, if and e are sufficiently small. Now, for 
n < K 2 1/4 (e + /3- 1 )- 1 /4 ) the 

sum is smaller than a constant multiplied by 

e + and this proves ([4"T[) . 

As for ([48 p . we observe that, since H is even in the momenta, F and H 
are uncorrelated, so that, using px,Y < 1, one gets 



|px„,if| < — 



— f l^-^i + ± \c Pj . H \] <^£^ + ?k^L 



As we have just shown, for n < n 2 (e+/9 1 ) the last term at the r.h.s. 
tends to zero as e + and in the same way behaves cr® 1 /ax n — 1- So, we 
limit ourselves to study ctq 1 -f/^q 1 = l/y^l + ^f/^Oi-F- By computing 
explicitly ct 2 ^ and applying the upper bound ([59[) to <r@ > <t© ^ we get 
that there exists a constant R > 1 such that 



^ei V k{s + (3- 1 ) 2 . 

Since the r.h.s. differs from 1 by a quantity larger than e 2 /3 2 , the corrections 
given by the other terms can be neglected if /3 > e^ 1 . This completes the 
proof. 
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Q.E.D. 



7 Relation between stability estimates and relax- 
ation times 

In the present section we discuss which implications the existence of an 
adiabatic invariant has in the frame of ergodic theory. The main point is 
that it can provide a lower bound to the relaxation time to equilibrium. 
Since there is no agreement in the literature on the definition of relaxation 
time, we will give here a mathematically clear form to such a concept. To 
this end, we need a preliminary discussion of an a priori bound to the time 
autocorrelations (see Theorem [3]). This is provided in Section [7.1[ Then, in 
Section 17.21 we define the concept of relaxation time. 

7.1 Relaxation times and time correlations 

One of the open problems in statistical mechanics is that of thermalization, 
i.e., to establish whether a system, starting from a given microscopic state, 
does attain thermodynamic equilibrium, and, if this is the case, to estimate 
the time scale needed to reach it. Such a time scale is usually called the 
relaxation time. From a physical point of view the situation is complicated, 
because certain systems, for example gases, reach equilibrium on a very 
short time scale, while others, for example glasses, are believed to reach 
equilibrium on geological time scales. 

Linear response theory (see pjJJ [20] ) shows that susceptibilities can be 
expressed in terms of the time autocorrelations of suitable dynamical vari- 
ables (namely, those conjugated to the perturbing field). In particular, the 
susceptibilities assume the equilibrium values only for measurements which 
last a time large enough, i.e., larger than the time needed by the time auto- 
correlations to become negligible. 

Now, the time correlations between pairs of dynamical variables are 
widely studied in the case of chaotic systems (see, for example, [211 G2] 
or the monograph [23]). For such systems, the correlations are known to 
tend to zero, as t — > oo, and one of the problems is to estimate the decay 
rate, for long times, of the time autocorrelations of all dynamical variables. 
This however amounts to give an upper bound to the time autocorrelations. 
From the standpoint of linear response theory, it is also significant to bound 
from below the time autocorrelations of suitably chosen dynamical variables, 
because this leads to a lower bound to the relaxation time. 
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Corollary [2] of Theorem [T] gives an estimate of such a kind, showing that, 
for the system here considered, the relaxation time is larger than a constant 
t, which is exponentially large in the perturbation parameters (and moreover 
does not depend on the number of degrees of freedom of the system). Results 
analogous to Corollary are of a general type. In fact one has 

Theorem 3 (Bound to the autocorrelation of a dynamical variable) 

Suppose that, for a dynamilcal variable X, there exists a constant rj > such 
that 

\\[X,H]\\< Wx ; (61) 

then one has 

C x (t)>l-\i?t 2 . (62) 

Remark. This theorem is a slight modification of Theorem 1 of |15j and 
is proved in the same way. On the other hand, we think that the decision 
to focus on the time autocorrelation, which we make here at variance with 
paper [15] , is crucial, if one aims at obtaining significant estimates in the 
thermodynamic limit. 

def 

Proof. Introduce the difference 5 = Xt — X . As Xt satisfies the Liouville 
equation and X is time-independent, one has dt5 = dtX t = —[H,Xt], which 
in terms of 5 takes the form 

d t S = -[H,S]+Y , (63) 

with Y = f —[H,X]. It is well known that, \i being invariant, the solutions 
of the Liouville equation are generated by a one-parameter group U(t) of 
unitary operators in the sense that X t = U(t)X. As 5(0) = 0, the solution 
of equation (l63j) is given by 

5= f U(t- s)Yds , 
Jo 

so that, U being unitary, one gets the estimate 



< 



o 



U(t-s)Y 



ds = t \\Y\\ < vjtax ■ 



Then, one gets the thesis by using the simple identity 

Cx(t) = 1 - . 



Q.E.D. 
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7.2 Definition and evaluation of relaxation times 

Taking into account the relation between susceptibilities and time autocor- 
relations, it is meaningful to introduce a parameter a = a(t), with values 
in [0, 1], which estimates how much the system is close to equilibrium, after 
a finite time t. So, for the set B C L 2 n C°° of the dynamical variables 
uncorrelated with H, we define 

Definition 1 (Correlation level) The correlation level a(t) at time t is 
defined as a(t) = sup^gg |Cx(i)|- 

Remark. One can limit oneself to the smooth observables, because these 
are the physically relevant onesll 

In the chaotic case a tends to as t — > oo: thus, looking for the decay to 
zero of the correlations is equivalent to looking at the asymptotic behaviour 
about zero of t(a), the inverse function of a(i)|3 On the other hand, accord- 
ing to linear response theory, it is more meaningful to look at the time after 
which the correlations are below a certain threshold. So, we introduce the 
following notion 

Definition 2 (Relaxation time relative to level a) The relaxation time 
relative to level a is defined as t(a) = inf t*(a), where t*(a) is such that 

sup \C x (t)\ < a for all t > t*(a) . 
XeB 

Remark. In order to provide a significant lower bound to the relaxation 
time t(a) at level a, as previously defined, it is clearly sufficient to find the 
time at which the autocorrelation of at least one dynamical variable X uncor- 
related with the Hamiltonian is certainly larger than a. Now, the following 
corollary on the relaxation time descends immediately from Theorem [3] 

Corollary 3 (Bound to the relaxation time) Suppose there exists a dy- 
namical variable X £ B and a constant rj > such that \\[X,H]\\ < rjax; 
then one has 

t(a) > y/2(l -a) - . 



Notice that one can find "pathological" functions for which the decay is arbitrarily 
slow (see (51]) even for strongly chaotic systems. For this reason, the control is usually 
restricted to a fixed continuity class. 

7 As a matter of fact, we cannot guarantee that a(t) is invertible, but we will give below 
a meaningful univocal definition of t(a) (see definition [2]). 
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The point is that many Hamiltonian systems of interest for Solid State 
Physics reduce to integrable ones in some limit, while, on the other hand, 
for integrable systems one has t{a) = +00 for any a < 1, since their integrals 
of motion remain correlated for all times. The question is then, what is the 
behaviour of such systems when the perturbation is small, i.e., to study the 
ergodic properties of slightly perturbed (or nearly integrable) Hamiltonian 
systems. It is natural to think that there exists a sort of continuity as the 
perturbation diminishes. Continuity can in fact be recovered in terms of the 
time needed for the system to reach thermalization (i.e., a sufficiently low 
correlation level). This is indeed the case in the system we have considered, 
because we can say that, as a consequence of Theorem HJ one has the lower 
bound 

t(a) > VW^a) HG^TF 1 )) ' 
which goes to infinity as both e and — > 0. 

8 Conclusions 

In this paper, we have constructed, for the Klein Gordon lattice, an adia- 
batic invariant, i.e., a dynamical variable whose time derivative is small as 
a stretched exponential with the perturbation parameters. Thus our result 
is similar to those which are known in Hamiltonian perturbation theory in 
the case of a finite number of degree of freedom or in the case of an infinite 
number of them, but at a fixed total energy (see \25\ 2]). The new feature 
of the present work is that our theorem remains valid in the thermodynamic 
limit, because the given bound turns out to be independent of the number 
of particles, and depends only on intensive quantities. As a corollary, we 
bound from below the stability time of such a model. 

We now add some comments. The first one concerns the fact that in our 
model we have two perturbation parameters, e and 1/(3. We believe however 
that the only really relevant parameter is Indeed, at least formally the 
parameter e can be arbitrarily decreased by performing a suitable normal 
form change of coordinates (|26J), and it seems to us that such a normal 
form does not alter in any fundamental feature the perturbation Hi (i.e., its 
local character). At the moment, however, we are unable to say anything 
definite on this point. 

In any case, while the estimate of Section [6] could presumably be applied 
to models more general than that studied in this paper, the estimates on 
the marginal probabilities of Section [5] are especially adapted to our model. 
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It would be important to improve our method, making it more flexible in 
order to cover more general situations. 

In our opinion, the real big problem that remains is the construction of 
adiabatic invariants for problems in which small denominators appear. This 
problem could be overcome in particular cases (see, for example, paper [IB]), 
but no precise strategy exists yet for the general case. We plan to tackle 
this problem soon. 
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A Estimates for the construction of the adiabatic 
invariant 

Here we intend to prove Lemma [T] of Section [H In order to do that, we need 
to recollect the usual algebraic properties used in perturbation theory (see 
[TT]). adapted to our norm || • || + , defined by (fTl"|h Such properties are stated 
in Lemmas [9l41 II later on, then the proof of Lemma [1] is briefly sketched. 

In order to develop the perturbation theory, a primary role is played by 
the action of the operator Lq and by the projections on its kernel and its 
range (see Section E|) , and these are more easily discussed in terms of the 
complex variables which diagonalize Lq. These are implicitly defined by 

qi = (6 + im) . pi = - im) > 1 < * < N (64) 

and in such variables one has 

L e V k =iu(\k\-\j\)er] k . (65) 

We must, however, take into account the fact that the norm || • || + is 
not invariant under such a change of coordinates. In fact, such a norm is 
formally well defined also for polynomials depending on the variables (£, rj) 
if, in the definition of H r s ' 1 and V s , r , we simply substitute for (p,q) the pair 
(^,rj). In that case, denoting by /' the transform of / via (|64p . one will 
have, in general, ||/||+ ^ ||/'|| + - On the other hand, the following lemma, 
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whose proof is identical to that of Lemma A.l of paper |17j . enables one to 
estimate the difference between the norms of the two functions. 

Lemma 9 Let f(q,p) be in V s ,r an d let /'((,, r/) be the transform of f via 
Then, one has f G V s ,r and ||/'|| + < 2l ||/|| + . Moreover, let g'(£,r,) 
be in V s r and let g(q,p) be the transform of g via the inverse of (64\). Then, 

' — l 

one has g G V s ,r and \\g\\ + < 22 \\g \\ . 

We need also the following lemmas 

Lemma 10 Let f be in V s ,r an d g in V s ' r '- Then, [f,g] G "P s + s '_2,r+r' an d 
one has, both in real and in complex variables, the inequality 

\\[f,g}\\ + < (2r + 2r' + l)ss' \\f\\ + \\g\\ + . 

Proof. See Lemma A. 2 of [T7j, noticing that, for any fixed i, each term 
of / contained in Til' 1 has Poisson bracket different from only with the 
monomials of 1-L r s , ,k such that \i— k\ < r + r'. The number of such monomials 
appearing in the decomposition of g is smaller than 2r + 2r' + 1. 

Q.E.D. 



Lemma 11 Let f G V s<r be a polynomial in complex variables. Then Hj\ff, 
Hjlf an d LQ 1 Hnf belong to V S r an d the following inequalities hold: 

||n^/|| + < ||/|| + , ||LW|| + < ||/|| + , ||Lo 1 iW|| + < 11/11+ • 

Proof. The fact that Lq 1 f belongs to V s ,r comes directly from Lemma [TOl 
as Hq is in Viq- The remaining statements are a consequence of the fact 
that Lq is diagonal in complex coordinates and that the smallest eigenvalue 
of Lq on 1Z has modulus oj > 1, in virtue of (l65|) . 

Q.E.D. 

Proof of Lemma[T]We pass to complex variables via Lemma[9]and proceed 
by induction on n, checking at each step even two supplementary inductive 
hypotheses: 

i) ^ n can be decomposed as ^ n = ^"=0 ■ > wriere *n £ ^21+2^-1; 
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ii) the following bound holds 



tt® < 2 n 2 10 ^ (n!) 2 (n - 1)1 , „ ,. | E ^' . 

+ l\(n — iy. 

As a matter of facts, on account of Lemma [TT1 such an estimate enables 
one to control the contributions due to Xn an d O n) which appear in the 
recurrent procedure that determines Xs, for s > n. Then, we come back to 
real variables via lemma [9] again. 

Q.E.D. 



B Technical proofs 
B.l Proof of Lemma [5] 

We start by proving formula (1341) . On account of the symmetry of the 
periodic system, one can pass from a system with N — l particles to one 
with N by inserting one more particle after the i—th site, for i = 1, . . . ,N— 1. 
The potential energy of the corresponding system is given by 



U N (qi, ...,qN,q) = Un-v 



2uj 



7i+i 



^2 , e 



B) +^ (q ~ qi + i) +^+ 



Neglecting the second term at the r.h.s. (which gives a contribution to the 
partition function which can be bounded from below by 1, and averaging 
over i in order to get a traslational invariant system, one gets 



Qn 

Qn-i -N 



N-l 



+oo 



8=1 
+ 0O 

dqexp 



dqi 

P 



dqN-i DN-i(qi, ■■■ , qN-i)x 



2uj 



7T~ W + j-+e{q 



2lu 



qi) 2 + e{q - q l+1 ) 2 



Here we have put qjy = q\. Then, we introduce the function <p qi (q) 
1 — exp \—/3e(q — qi) 2 /(2w)l , for which the inequality 



(66) 

dcf 



exp 



> 1 - ¥> 9i (g) - <Pq i+ i(q) 
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holds. We will show now that (f qi (q) is small except for a set of small 
measure. Making use of the previuos inequality, relation (|66p becomes 



Qn 
Qn- 



>a(f3,e) 



dqi . . . 



+00 



&q N -i £>N-i{qi, • • • , qN-i) x 



N-l 



/ 2 
~°° i=l 
in which the function a(/3,e) is defined bj|^| 



2u [ q + 2u ) 



(67) 



a{/3,e) 



def 



d</ exp 



'2w 



o 4 

! 2w 



2we» 



4 







where K a {x) is the Bessel modified function of second kind. The well 
known properties of K a (x) imply that a(f3,e) can be written as a(/3,e) = 
G((3, e)y / 27ru;//3, where G is a function always smaller than 1, approaching 
1, at fixed e, as j3 — )■ +00. We go on by dealing with the integral in ([67]) . 
first giving an upper bound for the innermost integral over q. We estimate 
it by splitting the phase space of the N — 1 particles periodic system in two 
sets: we will fix k > and consider Cl(N — 1, k), which is defined by 



Q(N-1, K ) = I (q 1 ,...,q N - 1 ) such that £ <g < ^-k(N - 1) \ , (68) 



N-l 



2lo 



i=l 



l3 



and its complement. In the latter set, the integral is simply bounded from 
above by 2a(/3,e). On the other hand, in order to estimate the integral in 
the set Cl(N — 1, k), we observe that, for any ki, the number of particles for 
which \qi\ > kik2u;//3 holds cannot exceed (N — 1)/k\. For these particles 
the integral is estimated again by 2a(/3,e). For the purpose of estimating 
the contribution of the remaining particles, we introduce the function 



I(/3,£,K,Ki 



def 



a(P,e) 



s up [ tp y (q)exjp (-^-qA dq 
/ k , k .9.,„/rJ-oc \ Zu) J 



\y\<y/^K2ufp 



We point out that /(/?, £. fay Kj\ ) tends to as e tends to 0, for (3, k, k\ fixed. 
Then, in the region Vt(N — 1,k), for any K\ > 1, one has the bound 



/+00 2 N ~ 1 
d ^T7— 7^2 exp 
-00 ^ 1 i=1 



l_ 

2oj 



2 , <f 
1 + 2^ 



< 



+ 2I(j3,£,K,Ki) 



a(j3,e) . 



Remark that the function a(/3, e) depends on e only via the term to = \/l + 2e. 



35 



We notice that we have provided estimates independent of qt , so the integrals 
over qi, . . . , qN-i appearing in (f67|) can simply be estimated as the product 
of these upper bounds times the measures of the sets in which the bounds 
hold. Now, we observe that the measure of £1 C (N — 1, k) is estimated by 



/ &q 1 ...dq N -ib N -i{qi,...,q N -i) < 

Jn c (N-i,K) 



Rn-i(P,k) 




(69) 



'f2 c (Af-l,K) WN-1 

where the function i?jy_i(/3, k) is defined by 

R N -i(/3, k) = f / dq 1 ... dq N ~! exp 

Jnc(N-i,K.) 

(Itt^^ fN-l , 

= (— ) r (— ■ K(JV_1) 

and r(s,x) is defined by (I42p . This way one obtains, finally, 

> ( 1 - - - 2J(A E , «, «x) - a(A e) . (70 ) 

Qn-i \ K i Qn-i J 



Prom this expression one can prove (JMJ) by induction on N. 

We now come to the proof of (|35p . We make use of the trivial inequality 
P(Ai n...ni r ) > 1 — Yll=i which holds for any probability and any 

collection of sets A\, . . . , A r . Consequently, we obtain 

Piv (l^l< A ... A \q r \<Q^j \ > 1-r-Pjv (\qi\>Q^y) ■ 

because, due to the translation invariance of the periodic system, every set 
has the same measure. Recall that Pjv^|<7i| > Qy/2u/P^ is just the integral 

of Dn times ^-| g .|>ey / 2^/^ " ^ ' 3 ° un< i to this integral can be found proceeding 
as above, i.e., by symmetrizing on qi, fixing k > and integrating separately 
over 0(iV, k) and its complement (recall that Q(N,k) is defined by (|68p ). 
This way we get 

(\qi\>&\f^f) < T^-Yl I 1 lnA>^/^7RD N {q 1 ,...,q N ) 

1 



Pj ~ NQ N ^J a{NtK) l«l>eVM* 
Rn(P,k) 



Qn 
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where Rjy is defined by ([69]) . and we bound -*-|, ? .|>e^/ 2aj//3 ^Y 1 m Q c (N,k). 
It is straightforward to notice that the number of sites for which > 
0a/2u;/ (3, in the interior of f2(iV, k), cannot exceed Nn/Q 2 . Therefore, the 
former term at the r.h.s. of the previous formula is smaller than l/4r if 
G > As far as the latter is concerned, we can choose k such that 

i?jv(/3, k)/Qn < l/4r, as we have shown above. For example, we can fix 
k = log(4rKo)- This suffices to infer that, for > 2y / rlog(4rKo), ([35]) is 
valid. 

B.2 Proof of Lemma [6] 

The first inequality in (j39[) comes directly from the fact that the integrand 
appearing in the definition of Qm is smaller than the function njv/,y, i-e., the 
integrand in the definition of Q\r- 

As regards the second inequality in ([39]) . we note that the integrand of 
Qm is equal to the one of Z M multiplied by p terms of the form exp(— f3eq mi q m 
at the sites, in number 2y, on the boundary of the blocks, which we denote 
by mi, . . . ,m-2 r , with the convention that m^y+i = m i- Then, we integrate 
only in the region in which the q coordinate of each of these sites is smaller 
than Qy/2u/p, with G = 2a/2j: log (8rK ), and we observe that 

Z M > exp (-AjceQ 2 ) Q M P M (\q mi \<Qy/2^ A ... A |g m2f | <e^JJ^j 

> ^ (8 ? Ko)- 32£o?2 ■ 

Here, Qm comes from the normalization of the probability, and in the second 
line use is made of Lemma 

We come now to inequalities (|40p and observe at once that the first one 
is trivial, because, on account of the identity in (i28j) . one has nM,f < n-M,p 
The second one is more complicated: we begin by proving it in the case in 
which each block is constituted by an even number of elements. 

In order to estimate Q M i we divide again the phase space of the system in 
the region £1 in which | q mi \ < y / 2un/(3, . . . , | q m2f \ < y2um/j3, where k > 
is a constant to be determined, and in its complement £l c . The integral over 
O is smaller than Qm ■ exp (4yeK) , while, as regards the complement, we 
notice that it is contained in the set in which Y^Li Qi — 2fiXc>//3. Thus, the 
integral over such a region is bounded from above by Q M -2% i^^/f^Y k), 
with fi < y, where we have dropped some positive term in the poten- 
tials, then we have integrated first over q mi , . . . , qm 2{ (which gives the term 
(2-KLoj ' P>fY{i, k), with r(y, k) defined by (j4*2j) ): then the blocks made of just 
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2 particles disappear, so that the integration over the remaining positions 
gives the term Q^_ 2 r This way we get 

(Into \ ' l 
-^J r(y,/c) • Q r A V_ 2y , with ti < y . 

Now, we apply the previous inequality to the function S^_ 2 r a t the 
r.h.s, and we end up with a relation similar to the previous one, in which 
however there appears the function QjW-2r-2ri' with r 2 < px- So, we can 
iterate this procedure, observing that r(y, k) is an increasing function of y, 
and we get 

Q l M < exp (4t£K) ^ f J (r(y, , with Q = f 1 , 

where we define Oj = Si=o ?fc> with yo = J, and J represents the integer such 
that a j = M/2. We make use of inequality (fM|) of Lemma [5] and finally get, 
if the series converges, Q f M < ew (^ £K ) Qm Yl'jLoiK^T^, n)y . We point 
out that the common ratio of this geometric series is a decreasing function 
of k, which tends to as n — > +oo: thus, we choose R = k(^,Kq) so as to 
satisfy (|4ip . and obtain the relation 

Q\f < 2 exp (AjceR{jc, K )) Q M ■ (71) 

If a number A < p of blocks is constituted by an odd number of elements, 
we integrate on one of the sites on the boundary of each of these blocks, in 
order that each of the blocks, in number r/, of the resulting lattice contains 
an even number of elements. By dropping some suitably chosen interaction 
terms in the potential, one gets Q\.j < (2vra;//3) A/2 Q r M _ A , with y' < y, where 
the blocks made of just one particles disappear. Now we can use (I7ip with 
Qm-\ instead of Qm- Then, making use of Lemma [5] to express Qm— a in 
terms of Q m , we get (l40|) . 

B.3 Proof of Theorem H 

As already said, the proof is performed by bounding from above every term 
at the r.h.s. of ([53]) . i.e., A(j,i), for j,i G V. 

We point out that the expectations in the definition of A(j, i) depend on 
the choice of Q, which is not completely fixed by its marginal probabilities 
(see comments on relations f|51 [) ) . In fact, the main part of paper [8] consists 
in introducing a suitable reconstruction operator (on the space of the joint 
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probabilities) which enables one to find a joint probability distribution that 
minimizes A, starting from an initially chosen one. We adopt the same 
technique, with the only difference that we apply it not at all sites of the 
lattice, but only at those lying on the complement of a fixed set V (we will 

call it f d = T\V). 

We also need to control, together with A(j,i), the auxiliary quantity 

where £ and £ 2 are the same Gibbsian fields entering in (|55|) . 

We introduce, then, the main tool of the proof, i.e., the reconstruction 
operator Ui, with i G T, which will enable us to construct the joint prob- 
abilities Q(dx, dy) on V (see formula ([54]) ). This operator is defined on a 
couple of fields (C 1 ,^ 2 ) having the same conditional probability at i, as fol- 
lows. For each pair of configurations x*,x 2 G M} T \, we denote by P\ x2 the 
measure on M 2 for which the minimum of the distance between P s x i and 
P; x 2 is attained, i.e., such that, for any measurable Bel, one has 

P xljx2 (M xB) = P i)Xl (S) , P xljx2 (P xR) = P i)X2 (P) , 
and / l x ^ y P^ lx2 (dx,dy) = D(P- i,P- 2) . 

Such a definition enables us to describe the action of Ui, because this op- 
erator maps the couple (^ 1 ,^ 2 ) into (C 1 ,^ 2 ) such that, for any measurable 
C C M 2 , 

P ((^>£?) e C I ^T\{i} = x T\{i}'^T\{i} = X T\{i}) = ^^i^) ' 

and, for any finite V C T not containing i, the joint probability distribution 
of (£,y,£,y) coincides with that of (Cy,^ 2 /). 

The effect of Ui on 7(i) and \(u, i) is described in detail in Lemmas 2, 3 
and 4 of the work [8]. Following such a paper, we adopt the convention that 
the quantities relative to the reconstructed couple are distinguished 

from the corresponding ones relative to (£ , £ 2 ) by adding the symbol For 
every set S C T we define the operator 

U s = U h oU i2 o...oU im , (72) 

where the order of the points ii, . . . , i m , contained in S U d^S is chosen in 
a suitable way. This is described in full detail in the proof of the following 



7(i) = E 
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Lemma H2J If we define 75 = sup ie s7(i) and As = supj i 6 s A(j, i), we can 
describe the action of [/^ on a couple of fields having the same conditional 
probability on SL)db r S accordingly to the following lemma, which is proved 
in Appendix IB, 41 

Lemma 12 Let (£ 1 ,£ 2 ) be a couple of fields having the same conditional 
probability on SUd^S, given by a specification T £ @(h, C, 6) flA(/t, KC, a), 
and (I 1 , | 2 ) = C/s(C\f 2 ). Then, one has 



where N, M and R are constants depending on a and b only. 

We remark that, if the eigenvalues of A are smaller than 1, the recon- 
structed quantities are smaller than the initial ones. So, we want to iterate 
the reconstruction procedure as much as possible. It turns out that we can 
iterate the procedure at most a number of times proportional to the distance 
between V and V. The reason is the following. 

In our case, £ is the field relative to the equilibrium Gibbs measure and 
£ 2 that relative to the probability conditioned to the configuration x on V, 
which we consider as fixed. It is apparent that such fields have the same 
conditional probability on every set which does not intersect V, but not on 
the whole T; by hypothesis, this conditional probability is that given by 
T G @(h, C, S) n A(/i, KC, a). Since the reconstuction procedure shrinks the 
set S on which we can control 7 and A, we can iterate it until V C S. So, 
the maximum number of iterations is attained if we start by reconstructing 
on V U d n brV, where n is the largest number such that d( n+ i^V D V = 0. 
We use Lemma [12] as the first step of a recurrent scheme, by applying each 
time Uv m , where V m+ i = V m \Jd^V m , Vq = V. In virtue of Lemma [T2| after 
the application of Uy m , one has 



Thus we get that the final values of 7^ and \y are smaller than the result 
of the application of the matrix A n to the vector with components 7, A. 




in which the matrix A is defined by 




a + NK- 1 C- l MK~ l 
C{R + NK- 1 ) 5a b + MK- 
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Moreover, we observe that we can write A = J 1 AJ, where A and J are 
defined by 

~. def / a + NK- 1 MR- 1 \ def ( C 

A " V i? + AT*" 1 «5a b + M^" 1 J and J " V 1 

This way we get A n = J~ 1 A n J. As the component Xy, which is the one we 
are interested in, is not affected by the action of J~ 1 , we can write that it 
is the second component of the matrix product 

X 

Since the eigenvalues of A are smaller than G *= max{(a + l)/2, (5a b + 
l)/2} < 1, if R is large enough, there exists Rq such that 

X v < (D/2) max{A, Cj}G n , 

where D is a constant depending on a, b, a, 5 only. On the other hand, 
n = d(V, V)/(br) and 7 < 1, from which there follows 

Ay < j{h)xexp(-cd(V,V) 



where c is defined in (|58p . 

In order to show (I57j) , we need only the use of (I54h in estimating the term 
in brackets of (|49p . We then observe that the r.h.s. of (|54p is smaller than 
2|F| 2 Ay, for the joint probability we have just found, and this concludes the 
proof. 



B.4 Proof of Lemma 1121 

Lemma 5 of work [8 J shows the result of the application of U{, in a suitably 
chosen order, to every site of T in sequence: one obtains that, for the couple 
of fields (£\£ 2 ), with the same specification T £ 0(/i, C, 5) n A.(h, KC, a) 
and R > 1, and for the reconstructed couple (^ 1 ,^ 2 ) the following matrix 
relation hold^l 

9 As a matter of fact, A is not the same matrix which appears in [8], since we needed 
to make the dependence on C explicit. Our statement can be proved by checking, in 
proving the induction (11)-(14) of [Sj, that the constants JV(-, •) are proportional to C 2 , the 
constants iV(-) and M(-, •) are proportional to C and the constants M(-) are independent 
of C. 
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In the proof of Lemma 5 of [8], the order of the f/j's is chosen in the 
following way: the lattice is partitioned in b disjoint sublattice, Zq, . . . , 
which are the cosets in T of 7L V with respect to Zq. Then the reconstruction 
is applied in sequence to each sublattice, and use is made of the fact that, 
if i € Zi, there exists a bound to 70), A(j,k) and A(k,j) for j G Zi\{i} and 
k G r\{i} which does not change after the application of U{. In particular, 
this implies that the bounds do not change for the already reconstructed 
sites in Z\. Neither does the reconstruction at site i change, on account of 
Lemma 4 of paper [8], the value of A(i,j) and A(j,i), for |j — i| > r. In this 
sense the reconstruction is local. 

So, the values of ^y and Ay, after one application of U\, depend at most 
on the values of 7(i) and A(j,i) in V U d r V. It is thus appearent that we 
can control the values of the reconstructed quantities only in a set V smaller 
than the set V on which we control 7 and A initially. In particular, V can 
be chosen so that V = V U d r V . Therefore, for any V C S, we define for 
I = 0, . . . , b— 1 a nested sequence of sets Vi + \ = f Vi U d r Vi, with Vq = V and 
the operator Uy, as Uy = Uy r,z • • • Uy b _ 1 f ] z b _ 1 , and notice that (see the 
above remarks) the order in which the sites in Vi (7 Z\ are chosen does not 
matter. Then, after the application of Uy, one has that 



for the same matrix A appearing in (173ft . This concludes the proof. 
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